Octahedral Tilt Instability of ReOs-type Crystals 
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The octahedron tilt transitions of ABX3 perovskite-structure materials lead to an anti-polar 
arrangement of dipoles, with the low temperature (T) structure having six sublattices polarized 
along various crystallographic directions. It is shown that an important mechanism driving the 
transition is long range dipole-dipole forces acting on both displacive and induced parts of the anion 
dipole. This acts in concert with short range repulsion, allowing a gain of electrostatic (Madelung) 
energy, both dipole-dipole and charge-charge, because the unit cell shrinks when the hard ionic 
spheres of the rigid octahedron tilt out of linear alignment. 



In 1950 Slater [1| presented an electrostatic theory of 
the ferroelectric transition in BaTiOa (perovskite ABX3 
structure) by generalizing the Clausius-Mossotti (CM) 
picture. An important ingredient is the fact that the 
local dipolar electric field Fjj at the X site differs from 
the Lorentz value (47r/3)P because the local symmetry is 
less than cubic. Here we present a similar discussion, also 
highlighting the role of dipole-dipole interactions, of the 
octahedron tilt commonly found in perovskite structure. 
To simplify the many interactions, this paper is confined 
to the ReOs structure type (3] with the perovskite A 
sublattice empty. We take AIF3 as our prototype. 

Starting with the parent cubic structure, the (i?3c 
rhombohedral) low T phase of AIF3 is generated by 
a rigid rotation of an AlFg octahedron through angle 
u! « O.Srad around a cubic (111) axis through an Al atom. 
Because neighboring octahedra share corners, rotations 
alternate, doubling the unit cell according to wavevec- 
tor (7r,7r,7r). A schematic view of the distorted (111) 
plane is shown in Fig. ^ The Shannon ionic radii 3] of 
AP+ and F~ match almost perfectly to an octahedron 
with the central Al touching the six F anions, and each 
F anion touching its eight F neighbors. Thus the rigidity 
of the octahedron follows both from Al-F covalency and 
from ionic size effects. As the octahedra tilt, their spacing 
shrinks, generating a rhombic primitive cell whose c/a ra- 
tio increases by the factor l/cos(a;) relative to the cubic 
value, and cell volume V decreases as AV/V — — sin^ uj. 

CM theory ^ j^] shows that a cubic lattice of polariz- 
able molecules may have a "polarization catastrophe" sig- 
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FIG. 1: The (111) planes of perovskite are alternately AX3 
and B layers. The AX3 sites constitute an fee lattice with 
close-packed triangular (111) planes occupied 75% by X an- 
ions (shown as filled circles) and 25% by A cations (shown 
as open circles and missing in the AIF3 structure.) Arrows 
denote displaced X anions. Faint triangles with counter- 
clockwise rotations have their nearest B cation in the plane 
below, and inverted faint triangles with clockwise rotations 
have their nearest B cation in the plane above. 



nailing an instability towards ferroelectric polarization. 
The condition for instability is an increase of the product 
na (density times polarizability) to 3/47r. A generalized 
version of this statement is derived in the Appendix: a 
self-stabilized spontaneously electrically polarized state 
will occur when a increases to l/7max where 7max is 
the maximum eigenvalue of the dipole-dipole interaction 
tensor T. The tensor Tiajp is defined as the a Carte- 
sian component of the dipolar electric field Fjj^i (at site 
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tilt angle CO (radians) 



FIG. 2: Electrostatic energy of BX3 with full ionic charges 
+3 and -1, as a function of the tilt angle u around the (111) 
axis. The dashed curve is the full quadratic approximation, 
namely the sum of the Madelung (M), dipole-dipole (D), and 
volume (V) terms. The solid curve is the exact Madelung 
sum. 

Ri of the lattice) created by a unit dipole in the Carte- 
sian [3 direction at site Rj. In a 3 — A^- vector notation, 
\F£i >— r|/i >. This generahzation of Clausius-Mossotti 
theory requires no restrictions on the size or symmetry 
of the system. The pattern of spontaneous polarization, 
given by the corresponding eigenvector, is in general not 
a simple ferroelectric. 

Rigid rotations of octahedra almost always occur in 
the low-T phases of perovskites 0, 01 • AIF3 has the per- 
ovskite structure with the A sublattice empty. Below 
730K the structure is rhombohedral because of cooper- 
ative rotation of AlFg octahedra, as shown in Fig. ^ 
We have successfully modelled this instability in two dif- 
ferent ways: (1) using density-functional theory (DFT) 
0, and (2) using classical molecular dynamics (MD) 9], 
including both short-range and electrostatic forces, plus 
anion polarizability. In this paper we abstract from our 
earlier results a simple picture that includes electrostatic 
effects of charged and polarizable point ions. 

The dipole patterns |/i > (< ia\fi > is the a Cartesian 
component of the dipole fli on the «'th ion) are particu- 
larly simple and illuminating. First define the "displacive 
dipole" flD,i of the ith X anion as jlD,i = —Peui where Ui 
is the small displacement from the site Ri of cubic sym- 
metry, and — /3e is the charge of the anion. The unknown 
dimensionless parameter /3 absorbs the uncertainty about 
what actual charge to assign. The dipole can be imagined 
as a separated pair of opposite charges, the negative end 
at the anion nuclear coordinate Ri + Ui, and the positive 
end at the lattice site Ri formerly occupied by the anion. 

Consider the dipole patterns generated by different 
types of octahedron tilts. In Glazer's f6| notation, when 



the tilt angle is small, any tilt (jjPip'>6'^ is a product of 
tilts around x, y, and z where p, q, r are -I- or -. There 
are six primitive rotations. For example, a z-tilt by an- 
gle 9 is denoted as cfPip^O^, where p ^ q ~ indicates 
no rotation around x or y. This tilt belongs to wavevec- 
tor k — (tTjTt, fc^). Rotations of adjacent octahedra in 
the X and y directions are forced to be opposite in sign 
(wavevector = ky = tt), whereas in the z-direction, 
the next octahedron can have a different rotation angle. 
Glaser's conventions are that r = + corresponds to kz = 
0, while r — — indicates k^ = n. We have discovered 
the interesting fact, explained in the Appendix, that the 
corresponding dipole patterns \fi > for the six primitive 
rotations are eigenvectors of F. The three \ii{(fPip'^6'^) > 
are degenerate with eigenvalue 7(-l-) = 14.383/a^, while 
the three \ij,{(j)'~'ip^9~) > are degenerate with eigenvalue 
7(-) = 14.461/a3 (where a = 3.43 A in AIF3). Note that 
7(— ) is the largest eigenvalue of F, and that 7(-l-) is only 
0.5% smaller than 7(— ). It follows that the dipole-dipole 
interaction energy of an arbitrary tilt-induced dipole pat- 
tern \n{(j)P^W) >=EC'7l7>is-E|C7p7/2- Thus an 
arbitrary Glazer tilt has dipole-dipole interaction energy 

Enir^en = -^(^) [7b)0'+7(9)V''+7M^'] 
= [7(pV'+ 7(9)^^+7 Wa^'] (1) 

where fj.^ = (}ea(j)/2 is the amplitude of the dipole eigen- 
array arising from the x rotations, etc. When p, q, r 
are all negative as in AIF3, the dipole-dipole interac- 
tion energy Eb = — 7£)|/2£)p/2 is as negative as possi- 
ble (7£) = 7(— ) is maximal) for any array of displacive 
dipoles of fixed magnitude For cases when the su- 

perscripts contain some -|-'s, i.e. for different Glazer tilt- 
ing schemes, the energy is at worst 0.5% smaller than 
optimum. 

So far we have presented an argument showing that 
tilts of corner-coupled rigid octahedra are strongly stabi- 
lized by dipole-dipole interactions. All such tilt systems 
are "6-sublattice" antiferroelectrics, because each of the 
six anions of a BXg octahedron defines a separate direc- 
tion of polarization. Although it is reasonable to char- 
acterize all such dipolar-stabilized tilt patterns as anti- 
ferroelectric, nevertheless conventional use often 0, ^ 
(but not always H'l) restricts "antiferroelectric" to cases 
where an applied external field can switch the state to 
ferroelectric. Therefore we use the term "anti-polar." 

Now we need a theory for the total energy. From our 
previous studies 0,0] it is clear that the displacive dipole 
/i D must be supplemented by an induced dipole /i/ on the 
polarizable anion. This costs energy ^^/2a per dipole, 
where a is the anion polarizability. The isolated anion 
in a field F has a moment aF which minimizes the to- 
tal energy ^'j/2a — fli ■ F. For the F~ anion, a depends 
somewhat on its environment, but is in the range [T^ 
^ 0.85± 0.05 A^. We also need to account for the to- 
tal change of electrostatic energy when anions move, not 
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just the cooperative dipole-dipole part. Taylor expanding 
the Madelung energy <i> to second order in displacements, 
two types of energy appear. (1) There is the intersite or 
dipole-dipole term already computed which alters the en- 
ergy bilinearly in the total moment /i/^^ + ^D,i of dipoles 
at different sites i and j. (2) There is a "Madelung" 
field Vi$ near Ri caused by the ideal undistorted lat- 
tice of other ions. This vanishes by inversion symme- 
try exactly at Ri and grows linearly {{ui ■ Vi)Vi$/2) 
with displacement from this site. This gives an energy 
^Ma/SlJ'DafJ'Dfs/'^- Uncertainty about actual charges is 
absorbed into the same factor /3 which was used in the 
definition of the displacive dipole fiD. From the Poisson 
equation V^<i> = 0, the trace of Tm is zero. For cubic 
BX3 with charges -1-3 for B and -1 for X, we find nu- 
merically that Tm is diagonal in cube body axes, with 
elements 7m = 40.789/0"^ in the directions transverse to 
the B-X-B axis (where rotations actually occur) and an 
element —2jm in the direction parallel. This describes a 
restoring field Fm = —^mUd at the site of the displaced 
X anion nucleus. This field is larger (by -2.8) than the 
destabilizing dipolar field = "fofJ-D previously found 
from the displacements of the other X anions. This ex- 
plains why the induced moment (/i/ — a{FM + Fjj)) is 
opposite to the displacive moments. 

The energy of the dipole array generated by a tilt is 

+ lAWDfJ.! - ^IVf^D- (2) 

The second, third, and fourth terms contain the Taylor 
expansion of the Madelung energy as described above. 
The fifth term accounts for destabilizing short-range in- 
teractions which have not yet been discussed. The large 
Madelung electrostatic attraction of ionic crystals tries 
to shrink the lattice constant as much as possible. Hard- 
core repulsion, a quantum effect, is needed to stabilize 
the lattice. In the simplified model of impenetrable hard 
spheres, the ions touch at the hard-sphere radius, and 
the lattice constant a of the cubic phase is twice the 
sum of the B and X ion radii. Pushing the X anions 
a distance u off the B-X-B axis by a tilt, the lattice 
constant changes to a -I- where Sa = —2u^ /a. The 
Madelung energy per cell of the B'^+Xj^^ cubic lattice 
is UMiV) = -17.908/32eVV"i/3 per ceh, where V ^ 
is the volume per cell. When V shifts to V + 5V , 
where 5V = a'^{Sax + Say + Sa^), the first order shift 
Um(V + SV) — Um{V) under a rigid tilt by u is thus 
-{l/2)jvfJ-D where jv = (4/3)f 7.908/0^ is the volume 
stabilization energy per dipole caused by the tilt u. 

The induced moment is the one which minimizes this 
energy, giving 

A^z/md = -(7M -7n)/(l/a-7D)- (3) 

For the parameters of AIF3, this gives fij ^ —O.SO^d- 
The total energy, evaluated at this optimal choice of in- 



duced moment, is 

C/tot(/^D, Ai/,opt) ^ 7^ [iM - ID -IV Tl ^J■D■ 

(4) 

If we had not included the anion polarizability (i.e. a — > 
0) then, for the parameters of B'^+Xg"^, the net restoring 
coefficient 7m — 7d — 7v would still be positive. In fact, 
our impenetrable sphere model probably overestimates 
7v', so stability is still fairly strong. However, there is 
a critical polarizability ac beyond which the quadratic 
restoring energy on dipoles goes negative, given by 

— =lD-\ ■ (5) 

ac IM - Id ~ iv 

For our simplified model of AIF3, this is 0.136A'^. Insta- 
bility thus occurs even for a well below the actual value 
for F-, a « 0.85 ± 0.05 

The theory presented here includes classical electro- 
static energies quite well (since non-overlapping charge 
distributions interact to good approximation as if they 
were point charges and dipoles) but lumps all quantum 
effects into a hard core, possibly overstating the amount 
of energy available in volume contraction. Our theory 
omits all higher than quadratic effects and thus cannot 
predict the magnitude of the tilt. Nevertheless, we have 
offered a sensible approximation with no free parameters, 
which helps explain nicely the nearly universal instabil- 
ity of perovskite materials to octahedron tilting. Fur- 
thermore, the model predicts very little energy discrimi- 
nation between different tilting schemes, consistent with 
the wide range of tilts seen experimentally (and some- 
times adopted by the same material at different temper- 
atures). The major influence of dipole-dipole interactions 
and anion polarizability are an interesting surprise. 
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APPENDIX 

A generalized CM theory can be constructed as 
follows. For an array of point polarizable molecules at 
fixed positions the energy for an arbitrary pattern of 
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FIG. 3: Eigenvalues of the dipole-dipole interaction F ver- 
sus wavevector for dipoles on the three F~ sublattices of the 
perovskite AIF3 structure. The largest eigenvalue is the five- 
fold degenerate state at the R point (tt, tt, tt), with eigenvalue 
14.461. At the right, the density of states is plotted horizon- 
tally versus energy vertical. 

dipoles {fli} is 

£ m}) = E - ^/^^ • ^^-) ^ ^ < Ml^i-r|/i > . 

(A.l) 

For dipoles on the X anions in perovskite structure, there 
are 3 sublattice sites per cubic unit cell. A vector space 
notation is used where |/i > is a 9A^-diniensional column 
vector of the 3 Cartesian components of each of the 3A^ 
dipoles in N unit cells, and T is the 9iV x 9A^ dipole-dipole 
interaction matrix which has elements 

3RijaRi]P - SapRfj 
i ia,j/3 = ■ y^-^> 

The lattice is stable against dipole formation if the 
quadratic form Ea. (|A.ip is positive (all its eigenvalues 
should be positive.) The condition for instability is that 



the maximum eigenvalue 7 of the matrix T exceeds the 
restoring force constant 1/a. The corresponding eigen- 
vector gives the pattern of displacement dipoles > 
which has the most self-stabilizing displacment pattern. 

Bloch's theorem allows eigenstates of F to be chosen 
as simultaneous eigenstates of translations T{R) where 
R is any translation vector of the primitive (simple cu- 
bic) lattice. The resulting Bloch states are labeled by 
wavevector k which lies in the first Brillouin zone. The 
eigenvalues were computed numerically, using an Ewald 
method to converge the sums. Results are shown in 
Fig. El As a test of the code, eigenvalues were also calcu- 
lated for a model with a fourth sublattice corresponding 
to the perovskite A sites. Together with the three X sites, 
the resulting lattice is equivalent to the face-centered cu- 
bic structure in a non-primitive conventional cube and 
a four-atom basis. It was found numerically that the 
largest eigenvalue equalled 47rn/3 (with n — 4/a'^) and 
occurred at fc = (000). This is the known result of CM 
theory - when every molecule bears the same moment pi 
and sits on a site of cubic symmetry, the field at each 
site is given by the classical Lorentz value (47r/3)P, and 
a ferroelectric instability occurs when the energy —P ■ F 
exceeds the cost fi^/2a to create the dipoles. 

The flatness of the uppermost branch of 7 versus k 
in Fig. El indicates that there is not much interaction 
coupling the xy-oriented dipoles in one xy plane to the 
xy-oriented dipoles in adjacent planes. This result can be 
understood as a consequence of the exponentially rapid 
transverse decay of electric field of a periodic array of 
dipoles This fact helps explain why the observed 

tilts of perovskites are so indiscriminating in their pre- 
ferred wavevector (7r,7r,0) vs. (tt, 7r,7r), and even allow 
mixed wavevector solutions. 
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